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Big Bang Nucleosyiitliesis represents perhaps the first, and still perhaps the most powerful partirle-astrophysirs 
connection. As such, it should provide an example for other work in this area. I discuss the current status of 
standard model BEN predictions and constraints, and then argue that the issue of observational systematic 
uncertainties is the key feature limiting our ability to constrain theory with observation. Nevertheless, several 
very important constraints are currently obtainable. For example, assuming maximal systematic uncertainties 
in ^He, D, and ^Li we find a conservative upper limit of fls < 0.16 in order for BBK predictions to agree with 
observations. Equally significant, we find that BBN predictions are inconsistent unless ''i/e abundance b)- mass 
is greater than 23.9%, or D -|-^He estimates are incorrect. By contrast, unless systematic uncertainties are taken 
into account, the quoted 2a observational upper limit on the primordial '^He fraction is 23.8%. 



1. THE POWER OF BBN 

Your honor, there are two facts in this case 
which are not in dispute: 

l.The predictions of standard BBN agree strik- 
ingly well with the inferred light element abun- 
dance estimates over many orders of magnitude 
2. As a result there exist strong constraints on 
both non-standard and standard particle physics 
and cosmology 

It is precisely these points which created the 
modern particle/astrophysics paradigm some 30 
years or so ago. We are now living in the 1990's 
however, and it is appropriate to move beyond 
self-congratulation over the successes of the stan- 
dard BBN model. In particular, there are two 
questions of particular relevance today: 

1. What exactly are the constraints? 

2. What precisely are the uncertainties? 

It is the second question which is of greatest in- 
terest here. I believe that it points at what will, or 
at least what should be the trends in this subject 
today. As cosmology turns more and more into 
an empirical science the issue of understanding 
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systematic uncertainties will become more and 
more central to utilizing data to constrain the- 
ory. Most important, if astroparticle constraints 
are to be believable, we must work hard not to 
over-interpret data, and derive constraints which 
may turn out later to be violated. As long as 
cosmology is dominated by observations and not 
experiments, the significance of systematic uncer- 
tainties should never be understated. 

Most of this review will focus on recent work 
performed at CWRU to address these issues[l-3]. 
Before proceeding. In the spirit of this meeting, 
I want to briefly remind the reader how powerful 
BBN constraints can be for particle theory. Not 
only can they constrain in principle the physics 
associated with processes at MeV-scale tempera- 
tures, they can allow us to have a lever- arm which 
could in principle extend to the Planck scale! A 
few examples should suffice: 

1.1. The Gramtino Problem 

This problem has been around in one way or 
another as long as local supersymmetry has. The 
gravitino is the spin 3/2 supersymmetric part- 
ner of the graviton. In all such models of cur- 
rent interest, supersymmetry breaking is com- 
municated to the observed particles by gravity, 
and the scale of supersymmetry breaking is man- 
ifested by a gravitino mass. Since supersymmetry 
fixes the gravitino couplings to matter to be iden- 



2 



tical to the graviton couplings, the gravitino can 
decay into all ordinary particles with a coupling 
strength which is purely gravitational, and hence 
very weak. As a result, the graviton decay rate 
is fixed to be of order F « (M^/M^i) where M is 
the gravitino mass. This implies that the gravi- 
ton lifetime is generally much longer than 0(1 
sec) unless M > (OlOTeV). Because the grav- 
itino is coupled to all particles, its decays will 
produce particles which can in turn photofission 
deuterium and helium, and thus destroy the good 
agreement of BBN predictions with observations. 
This problem can be resolved by diluting primor- 
dial gravitinos by inflation, but stringent require- 
ments on reheating scales result (eg. [4,5]). 

1.2. Another Technicolor Problem 

While technicolor theories arc beset with prob- 
lems, they do have the virtue of being the only 
theories to attempt to fully dynamically address 
the problems of the observed quark and lepton 
mass matrices. If technicolor proports to be a the- 
ory of mass, it must also address, at some point, 
the possilnlity of neutrino masses. This is nor- 
mally a difficulty for technicolor models, because 
if a see-saw type mechanism is utilized to pro- 
duce Majorana masses, the smallness of the tech- 
nicolor, or extended technicolor mass scale im- 
plies that the light neutrinos should not be very 
light. One way aroimd this is to postulate light 
Dirac masses, which can be achieved by postulat- 
ing symmetries which protect the observed neu- 
trinos from heavy Dirac mass terms which give 
masses to techni-neutrinos. However, perhaps the 
most stringent constraint on such model building 
comes from BBN. The extra-right handed neu- 
trinos are effectively sterile. Nevertheless, their 
interactions with their left-handed partners, re- 
sulting from extended tocliicolor interactions, are 
sufficient to populate them in the early universe 
at an unacceptable level unless the ETC scale is in 
excess of a few TeV [27]. This severely constrains 
such models. 

1.3. Neutrinos 

BBN is famous for constraining neutrinos. As I 
shall describe, this situation is getting even more 
interesting. We are on the verge of being able to 



address whether even completely "sterile" light 
neutrinos are allowed by BBN. In a related vein, 
even a singlet Majoran may be in trouble. 

1.4. Dark Matter 

BBN provides the strongest evidence avail- 
able that the universe is not closed by Baryons. 
In turn, BBN constraints also provide strong 
motivation for the possibility that the inferred 
dark matter dominating galactic dynamics is non- 
baryonic. As I shall describe, results allow the 
possibility in principle to convincingly resolve 
this issue. At the same time, recent results on 
the baryon fraction of dense clusters provides a 
potential confrontation with the idea that the 
universe is flat, unless systematic uncertainties 
forced BBN constraints to be relaxed. With this 
in mind, several groups have recently investigated 
how large the baryon abundance of the universe 
may be pushed in standard BBN models. 

2. A BRIEF PRIMER 

Since we shall be interested here in how BBN 
predictions depend on several fundamental micro- 
physical and cosmological parameters, including 
the number of light neutrinos, the neutron life- 
time, and the baryon density of the universe, I 
thought I would spend a few paragraphs outlin- 
ing why this is the case. 

First, the heart of BBN calculations is the 
equation which is at the heart of much of particle- 
astrophysics, the Boltzmann equation in an ex- 
panding universe. The number density, n, of any 
particle species evolves as 

dn/dt = - < crv >a n^ + Si< (jv >^ Xf-3nR/R 

Here, the first term on the RHS represents the 
thermally averaged annihilation cross section, the 
second term represents the sum over thermally 
averaged production cross sections, and the final 
term represents dilution due to the expansion of 
the universe. 

This equation permeates BBN considerations. 
It is responsible for (a) determining the remnant 
neutron number density, (b) determining the rem- 
nant abundances of all light elements, and (c) 
determining the number density of light neutri- 
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nos and other exotic particle species which might 
affect the expansion rate of the universe during 
BBN. The key point is tliat once tlio last term in 
the equation dominates the two preceeding terms, 
the number density of a given particle species be- 
gins to depart from its thermal orjuilibrium value. 
Massless particles continue to keep their thermal 
equilibrium form, but the temperature character- 
ising the distribution need not be the same as the 
temperature of the remaining particles in thermal 
equilibrium if entropy subsequently gets dumped 
into the radiation gas (Hence, some particles can 
have an abundance which is equivalent to a frac- 
tion of a neutrino in thermal equilibrium). Mas- 
sive particles, such as neutrons, have a distribu- 
tion which freezes out, so that their number den- 
sity is suppressed by a Boltzmann factor charac- 
teristic of the temperature at freeze-out, rather 
than continuing to fall with decreasing temper- 
ature. Finally, composite particles such as he- 
lium nuclei will not be produced with their ther- 
mal equilibrium value until the abundance of re- 
actants involved in their production (i.e. Xi) 
reaches a critical value. 

This completely explains the dependence of the 
residual ^He abundance produced during BBN on 
the fundamental parameters described briefly in 
the subsections below. Since one of the purposes 
of this review will be to demonstrate that almost 
all the BBN action today revolves around ^iJe, 
this will suffice for our needs. 

2.1. The neutron lifetime 

The neutron lifetime is one way of parametriz- 
ing the strength of the weak interactions which 
interconvert neutrons and protons. The longer 
the neutron lifetime, the weaker these reactions 
are. As a result, the longer the neutrino lifetime, 
the earlier the weak interactions which keep neu- 
trons and protons in thermal equilibrium decou- 
ples. As a result, more remnant neutrons will be 
available to partake in BBN reactions. In turn, 
more helium can be produced during BBN. 

2.2. The number of neutrinos 

The expansion rate of the universe is directly 
proportional to the density of the radiation gas, 
which is directly proportional to the number of 



species in the radiation gas. As this number in- 
creases by an amount which is equivalent to one 
extra neutrino helicity degree of freedom, the ex- 
pansion rate increases by a fixed amount. This in 
turn implies that the weak interactions decouple 
at a slightly higher temperature, which in turn 
results in more primordial ^iJe being produced. 
One of the recent results I shall describe is a new 
derivation of the relation between primordial he- 
lium and neutrino number. 

2.3. The haryon density 

Primordial ^He cannot form until enough deu- 
terium forms so that helium production reactions 
can compete with the expansion rate of the uni- 
verse. The earlier this happens, the more efficient 
is helium production (in part because neutrons 
are decaying slowly until they get bound in nu- 
clei). The greater the primordial baryon density, 
the greater the density of neutrons and protons. 
The greater the density of neutrons and protons, 
the larger the production rate of deuterium. The 
larger the production rate of deuterium, the ear- 
lier a critical density of deuterium nuclei forms 
which can result in rapid production of helium. 
Thus, the greater the baryon density the larger 
the primordial helium fraction produced. 

3. BBN: RECENT DEVELOPMENTS 

The remarkable agreement of the predicted pri- 
mordial light element abundances and those in- 
ferred from present observations yields some of 
the strongest evidence in favor of a homogeneous 
FRW Big Bang cosmology. Because of this, sig- 
nificant efforts have taken place over 20 years to 
refine BBN predictions, and the related obser- 
vational constraints. Several factors have con- 
tributed to the maturing of this field, including 
the incorporation of elements beyond ''He in com- 
parison between theory and observation [6], and 
more recently: an updated BBN code [7]. a more 
accurate measured neutron half life[8], new esti- 
mates of the actual primordial ''He , D +^ He, 
and 'Li abundances [9.10]. and finally the deter- 
mination of BBN uncertainties via Monte Carlo 
analysis [11]. All of these, when combined to- 
gether[12], yield a consistent and strongly con- 
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strained picture of homogeneous BBN. 

We recently returned to re-analyze BBN con- 
straints initially motivated by three factors: new 
measurements of several BBN reactions, the de- 
velopment of an improved BBN code, and finally 
the realization that a correct statistical determi- 
nation of BBN predictions should include correla- 
tions between the different elemental abundances. 
Each serves to further restrict the allowed range 
of the relevant cosmological observables Qb and 
Np. Of course, statistically determined uncer- 
tainties are not the major factor limiting our abil- 
ity to use BBN to constrain fundamental param- 
eters. As we shall see, systematic uncertainties in 
the inferred light element abundances are gener- 
ally much larger, and must be properly accounted 
for if we are to conservatively compare predictions 
with observations. In this first section I outline 
the details of our effort to properly update and ac- 
count for BBN statistical uncertainties, and leave 
the discussion of systematics to a later section. 

3.1. New BBN Reaction Rates: 

By far the most accurately measured BBN in- 
put parameter is the neutron half-life, which gov- 
erns the strength of the weak interaction which 
interconverts neutrons and protons. Since this 
effectively determines the abundance of free neu- 
trons at the onset of BBN, it is crucial in de- 
termining the remnant abundance of "^He . With 
the advent of neutron trapping, the uncertainty in 
the neutron half-life quickly dropped to less than 
0.5/^ by 1990. Nevertheless, it is the uncertainty 
in this parameter that governs the uncertainty in 
the predicted ''He abundance. The world average 
for the neutron half-life is now r,v = 887 ± 'Iscc 
[8] (Note: This is a recent update from the ear- 
lier value of Tjv = 889 ± 2.1 sec. Unless otherwise 
stated, the new value is used in the tables and for- 
mulae presented here and in our most recent work 
[3], which thus updates values found in some of 
our earlier work[1.2]). This has an uncertainty al- 
most twice as small as that used in vaious earlier 
analyses. 

We also updated the rate ' Be + p ^ +^ B 
[13] 20% which we thought might be signifi- 
cant at high values of rjiQ (defined by Qb = 
.00367i-2(T/2.726)^?7io x 10", where T is the 



microwave background temperature today, and 
h defines the Hubble parameter H = lOO/i 
km/(Mpc sec)), but is not. Other than these two 
new rates we used the reaction rates and uncer- 
tainties from [12]. 

3.2. New BBN Monte Carlo: 

Because of the new importance of small cor- 
rections to the ''He abundance when comparing 
BBN predictions and observations, increased at- 
tention has been paid recently to effects which 
may alter this abundance at the 1% level or less. 
In the BBN code several such effects were incor- 
porated, resulting in an rjiQ -independent correc- 
tion of -I-.0006 to the lowest order value of 1^ (the 
^He mass fraction). This is a change of -I-.0031 
compared to the value used in previous published 
analyses[9,ll]. 

In the present code, more than half of the new 
correction is due to finer integration of the nuclear 
abundances. Making the time-step in the code 
short enough that different Runge-Kutta drivers 
result in the same number for the ''He abundance 
produces a nearly ijiq independent change in Yp 
of +.0017 [15]. Residual numerical uncertainties 
are small [15,16]. The other major change is the 
inclusion of A/y^ effects[17]. Seckel showed that 
the effects on the weak rates due to nucleon recoil, 
weak magnetism, thermal motion of the nucleon 
target and time dilation of the neutron lifetime 
combine to increase Yp by ~ .0012. Also included 
in the correction is a small increase of .0002 in Yp 
from momentum dependent neutrino decoupling 
[18,19]. 

Finally, we utilized a Monte Carlo procedure 
in order to incorporate existing uncertainties 
and determine confidence limits on parameters. 
Such a procedure was first carried out [11] with 
BBN reaction rates chosen from a (temperature- 
independent) distribution based on then existing 
experimental uncertainties. This procedure was 
further refined[12] by updating experimental un- 
certainties and using temperature dependent un- 
certainties. Here we utilized the nuclear reaction 
rate uncertainties in [12] (including the temper- 
ature dependent uncertainties for ^He(a, 7)^Be 
and ^H(a,7)^Li) except for the reactions we up- 
dated. Each reaction rate was determined using a 



5 



Gaussian distributed random variable centered on 
unity, with a 1 — <t width based on that quoted in 
[12]. For tlio rates witliout temperature- depen- 
dent uncertainties this number was used as a mul- 
tiplier throughout the integration. For the two 
rates with temporaturo-dopondont uncertainties 
the original uniformly distributed random num- 
ber was saved and mapped into a new gaussian 
distribution with the appropriate width at each 
time step. 
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Figure 1. BBN Monte Carlo predictions as a func- 
tion of r\\Q. Shown are symmetric 95% confidence 
limits on each elemental abundance. Also shown 
are claimed upper limits inferred from observa- 
tion. 



The results of the updated BBN Monte Carlo 
analysis are displayed in figure 1, where the sym- 
metric 95% confidence level predictions for each 
elemental abundance are plotted. Also shown are 
previously claimed observational upper limits for 
each of the light elements [9.10.12]. In tlie first 
instance we shall utilize these limits in order to 
assess how BBN constraints have evolved based on 
our re-analysis, and after this we shall consider 



how to account for systematic uncertainties. Fig- 
ure 1 also allows one to assess the significance of 
the corrections we used, in relation to the width of 
the 95% C.L. band for Yp, which turns out to be <^ 
.002. The total change in Yp of « -I-.003 from pre- 
vious BBN analyses conspires with the reduced 
uncertainty in the neutron lifetime, which nar- 
rows the uncertainty in Yp and feeds into the un- 
certainties in the other light elements, to reduce 
the range where the predicted BBN abundances 
are consistent with the inferred primordial abun- 
dances. 

3.3. Statistical Correlations Between Pre- 
dicted Abundances: 

While the introduction of a Monte Carlo proce- 
dure was an important step, the determination of 
limits on the allowed range of BBN parameters 
Q,B and based on comparison of symmetric 
95% confidence limits for single elemental abun- 
dances with observations, as has become the stan- 
dard procedure, overestimates the allowed range. 
This is because the BBN reaction network ties to- 
gether all reactions, so that the predicted elemen- 
tal abundances are not statistically independent. 
In addition, the use of symmetric confidence lim- 
its is too conservative. Addressing these factors 
is a central feature of our work. 

Figure 2 displays the locus of predicted val- 
ues for the fractions Yp and D -|-^He/H for 1000 
BBN models generated from the distributions 
described above for rjiQ = 2.71 (figiu'e a) and 
rjiQ = 3.08 (figure b). Also shown is the = 4 
joint confidence level contour derived from this 
distribution, in a Gaussian approximation, calcu- 
lating variances and covariances in the standard 
manner. The horizontal and vertical tangents to 
this contour correspond to the individual sym- 
metric ±2(7 limits on Gaussianly distributed x 
and y variables. As is evident from the figure, 
and as is also well known on the basis of analyt- 
ical arguments, there is a strong anti-correlation 
between Yp and the remnant D -|-^He abundance 
(the normalized covariance ranges from -0.7 to - 
0.4 in the rjiQ range of interest). Thus, those mod- 
els where ''He is lower than the mean, and which 
therefore may be allowed by an upper bound of 
24% on Yp, will also generally produce a larger 
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remnant D+^He/H abundance, which can be in 
conflict with the bound on this combination of 
10~^ [24]. This will have the effect of reducing 
the parameter space which is consistent with both 
limits. 
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Figure 2. Monte Carlo BBN predictions for Yp vs 
D +^ He and allowed range for (a) r]io = 2.71, 
and (b) rjiQ = 3.08 (using 1992 neutron half life 
value). In (a) a Gaussian contour with ±2(7 limits 
on each individual variable is also shown. 



Because our Monte Carlo generates the actual 
distribution of abundances. Gaussian or not, we 
determined a 95% confidence limit on the allowed 
range of ijioiN,^) by requiring that at least 50 
models out of 1000 lie within the joint range 



bounded by both the ^He and D +^He upper 
limits, as shown in figure 2. This is to be com- 
pared with the procedure which one would fol- 
low without considering joint probability distri- 
butions. In this case, one would simply check 
whether 50 models lie either to the left of the D 
+^He constraint for low rjio (figure a), or below 
the ^He constraint for high rjio (figure b). This is 
clearly a looser constraint than that obtained us- 
ing the joint distribution. Finally, the procedure 
which has been used to-date, which is to check 
whether the symmetric 2(7 confidence limit (i.e. 
when 25 models exceed either bound) for a single 
elemental abundance crosses into the allowed re- 
gion gives even a looser constraint, as can be seen 
in figure 2a. 

In table 1 we display these results. Here we 
show the 95% confidence limits on ij^q, both as 
we have defined them above and also using the 
looser procedures which ignore correlations. Note 
that the constraints tighten dramatically as the 
number of effective light neutrino species, is 
increased. 

3.4. Implications and Caveats: 

The above constraints on rjio and N^, taken 
at face value, assuming a Yp upper bound of 24% 
and an D-|-'^He/H upper limit of 10^"*, would have 
significant implications for cosmology, dark mat- 
ter and particle physics. The limit on rjio corre- 
sponds to the limit 0.015 < Or < 0.070. (As- 
suming 0.4 < h < 0.8, as is required by direct 
measurements and limits on the age of the uni- 
verse.) Thus, if the previously quoted observa- 
tional upper limits on the Yp and D+^He/H are 
used directly, homogeneous BBN would imply: 

(a) The upper limit on fls would be marginally 
incompatible with even the value of « 0.1 inferred 
from rotation curves of individual galaxies, fur- 
ther suggesting the need for non-baryonic dark 
matter in these systems. 

(b) The bound on the number of effective light 
degrees of freedom during nucleosynthesis is very 
severe, corresponding to less than 0.04 extra light 
neutrinos with the 1992 neutron half life and .07 
extra neutrinos with the new world average (see 
later section). This is perhaps the most worri- 
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some bound of all because it argues against any 
Dirac mass for a neutrino without some signif- 
icant extension of the standard model. This is 
because even a light "sterile" right handed com- 
ponent whose interactions freeze out about 300 
GeV will still contribute in excess of 0.047 ex- 
tra neutrinos during BBN [20] without extra par- 
ticles whose annihilations can further suppress 
its abundance compared to its original thermal 
abundance. It can have further implications for 
a variety of kinds of hot or cold dark matter. For 
example, new light scalars would be ruled out un- 
less they decouple above the electroweak scale, as 
would be any significant population of supersym- 
metric particles during BBN. Even allowing 0.047 
extra light neutrinos, the upper limit on a Dirac 
mass would be reduced to « 5 keV [25,26]. A v^- 
mass greater than 0.5 MeV with lifetime exceed- 
ing 1 sec. would also be ruled out due to its effect 
on the expansion rate during BBN[21,23]. Also, 
neutrino interactions induced by extended tech- 
nicolor at scales less than O(IOO) TeV are ruled 
out [27], and sterile right handed neutrinos [28] 
would be ruled out as warm dark matter as the 
lower limit on their mass would now be O(lkcV). 
(c) The primordial '^He mass fraction would have 
to he greater than 23.9 % for conststency with 
D+^He constraints! [3] As we shall describe, 
this generally exceeds the 2cr upper limit on the 
primordial ''He abundance inferred without ac- 
counting for possible systematic uncertainties in 
the analysis. 

These constraints are so stringent that they cry 
out for a consideration of uncertainties in the light 
clement abundance estimates. Lidoed. as we shall 
next discuss, in spite of the considerable effort 
devised above to accomodate statistical uncer- 
tainties in the predictions, by far the largest and 



most significant uncertainties in the comparison 
of BBN predictions and observations come from 
the latter. Moreover, these uncertainties are sys- 
tematic and not statistical. Accomodating them 
in a BBN analysis will be the subject of the rest 
of this review. 

4. ABUNDANCE ESTIMATES AND 
SYSTEMATIC UNCERTAINTIES 

As we have just indicated, the weakest link in 
a BBN analysis involves the assumed light ele- 
ment abundances. Estimates of ^He , for example 
[30-32] are mostly indirect, and subject to large 
systematic uncertainties, which may also be im- 
portant for the other abundance estimates. As a 
result, our refined BBN analysis described above 
suggests the need for revision of the light element 
abundance estimates inferred from observation at 
least as much as it argues for or against new non- 
standard physics. 

One of the most worrisome aspects of the 
present constraints is the fact that D+^He pro- 
vides a lower limit on Y'p which is uncomfortably 
close the previously claimed upper limit. There 
are obviously two ways out of this dilemma: ei- 
ther the D+^He limit is increased, which would 
allow smaller values of Yp to be consistent, or the 
observational upper limit on Yp is increased. Both 
possibilities have recently been discussed as I shall 
describe below. In addition, recently, several 
groups have assessed more carefully the system- 
atic uncertainties present particularly in the pri- 
mordial ^He abundance estimates[36,37,42], and 
have fjuoted various new upper limits on cosmo- 
logical parameters based on their assessments. It 
is very clear, based in part on the differing esti- 
mates, that it is quite difficult at the present time 
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to get an accurate handle on these uncertainties. 

Because of this, and because we could utilize 
the full statistical machinery we previously de- 
veloped when comparing predictions to "obser- 
vations", we felt it would be useful to prepare 
a comprohonsivo table of constraints on N,^ and 
for a relatively complete range of different as- 
sumptions about light element abundances. This 
allowed us to explore the role of different esti- 
mates in the constraints, as well as the effect of 
correlations as the light element abundance esti- 
mates vary. In addition, it allows us to address 
several points which we feel are important to con- 
sider when deriving cosmological constraints us- 
ing BBN predictions. Finally, this analysis leads 
to new simple relations between the light element 
abundances and limits on cosmological parame- 
ters such as the number of neutrinos , N^, and 
the baryon to photon ratio rjiQ. 

Before proceeding, one might wonder whether 
if systematic uncertainties in the inferred primor- 
dial element abundances are dominant, one need 
concern oneself with the proper handling of statis- 
tics in the predicted range. Thankfully, there are 
two factors which make the comparison of predic- 
tions and observations less ambiguous in the case 
of BBN: 

(1) Because the allowed range in the observa- 
tionally inferred abundances is much larger than 
the uncertainty in the predicted abundances, any 
constraint one deduces by comparing the two de- 
pends merely on the upper or lower observa- 
tional limit for each individual element, and not 
only both at the same time. Thus, one is not so 
much interested in the entire distribution of al- 
lowed abundances as one is in one extremum of 
this distribution. 

(2) Systematic uncertainties dominate for the 
observations, while statistical uncertainties dom- 
inate for the predictions. 

Both of these factors suggest that a conserva- 
tive but still well defined approach involves set- 
ting strict upper limits on Yp, D+^He, and ^Li, 
and a lower limit on D, which incorporate the 
widest range of reasonably accepted systematic 
uncertainties. Determining what is reasonable in 
this sense is of course where most of the "art" 
lies. Nevertheless, once such limits are set and 



treated as strict bounds, then one can compare 
correlated predictions with these limits in a well 
defined way. In this way one replaces the ambi- 
guity of properly treating the distribution of ob- 
servational estimates with what in the worst case 
may be a somewhat arbitrary determination of 
the extreme allowed observational values. 

Clearly all the power, or lack thereof, in this 
procedure lies in the judicious choice of observa- 
tional upper or lower limits. Because of our con- 
cern about the ability at present to prescribe such 
limits I consider a variety of possibilities here. 
Once one does choose such a set, however, it is 
inconsistent not to use all of it throughout in de- 
riving constraints. If one uses one oliservational 
upper limit for Yp, for example, to derive con- 
straints on the number of neutrinos, but does not 
use it when deriving bounds on the baryon den- 
sity, then probably one has not chosen a suffi- 
ciently conservative bound on Yp in the former 
analysis. It has been argued that a weak, loga- 
rithmic, dependence of Yp on rjiQ invalidates its 
use in deriving bounds on the latter quantity. Not 
only can this argument be somewhat misrepresen- 
tative for an interesting range of Yp values, but 
until Yp exceeds statistically derived upper lim- 
its by a large amount, it can continue to play a 
signifcant role in bounding rjio from above. 

4.1. Abundance Estimate Uncertainties: 
The Range 

(a) ^He: By correlating ''He abundances with 
metallicity for various heavy elements including 
0,N and C, in low-met allicity HII regions one can 
attempt to derive a "primordial" abundance de- 
fined as the intercept for zero metallicity. This 
can be determined by a best fit technique, as- 
suming some linear or quadratic correlation be- 
tween elemental abundances (i.e. see [22,44-46]). 
The statistical errors associated with such fits 
are now small. Best fit values obtained typically 
range from .228-. 232, with statistical "Irr" errors 
on the order of .003-. 005. This argument yields 
the upper limit of .24 [22] which has been oft 
fluoted in the literature. Eecently this number 
has begun to drift upwards slightly. New ob- 
servations of HII regions in metal poor galaxies 
have tended to increase the statistically derived 
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zero intercept value of Yp by perhaps .005 (i.e. 
[36,47]). In addition, the recognition that system- 
atic, and not statistical uncertainties may dom- 
inate any such fit has become more widespread 
recently. The key systematic uncertainty which 
interferes with this procedure is the uncertainty in 
the ^He abundance determined for each individ- 
ual system, based on uncertainties in modelling 
HII regions, ionization, etc used to translate ob- 
served line strengths into mass fractions. Many 
observational factors come into play here (sec [48] 
for a discussion of observational uncertainties), 
and people have argued that one should add an 
extra systematic uncertainty of anywhere from 
.005-. 015 to the above estimate. Clearly thus, 
one should examine implications of ^He abun- 
dances in the range .24-. 25. We shall show that 
for Yp above .25: (a) ''Ho becomes unimportant 
for bounding 7710, and (b) the effect on bounds on 
can be obtained by straightforward extrapola- 
tion from the data obtained for the range .24-. 25. 

(b) ^Li: It is by now generally accepted that 
the primordial abundance of ^Li is closer to the 
Spite Pop II plateau than the Pop I plateau. Nev- 
ertheless, even if one attempts to fit the primor- 
dial abundance by fitting evolutionary models to 
the Pop II data points[10]. assuming no deple- 
tion, one still finds a 2(7 upper limit as large as 
2.3 X 10^^". The role of rotationally induced de- 
pletion is still controversial. It is clear some such 
depletion is expected, and can be allowed for[41], 
but observations of ^Li, which is more easily de- 
pleted, put limits on the amount of ^Li depletion 
which can be allowed. We assumed an extreme 
factor of 2 depletion as allowable, and thus we 
explore how cosmological bounds are affected by 
a 'Li upper limit as large as « 5 x 10^^". 

(c) D and D+^He: The situation regarding 
this combination has become quite interesting re- 
cently. There has been a new claimed observa- 
tion [33], of deuterium in a primordial gas cloud, 
at a level (D-F-^He)/H = 1.9 - 2.5 x lO"'*. It 
has long been argued that any present measure- 
ment of D provides a lower limit on its primor- 
dial abundance because D is so fragile that it is 
easily destroyed in stars. Also, because the pre- 
dicted BBN abundance falls monotonically with 
increasing baryon density, a lower limit on deu- 



terium thus places a reliable upper limit on the 
baryon density of the universe. Previously quoted 
solar system abundance estimates of 10"'^ led to 
a firm upper bound on 7710 < 8 which clearly es- 
tablished that baryons could not close the uni- 
verse. The Songaila ct al. observation, an or- 
der of magnitude larger, is also a factor of two 
greater than the previous upper limit on the com- 
bination D-f'^He. As a result, this would allow 
smaller values of jjiq, which would in turn allow 
a smaller value of Yp, although the upper limit 
on Q,B one might derive would be much more se- 
vere. We have recently explored the implications 
of this possibility in some depth [2]. In particular, 
we have shown that this result, if upheld, would 
change the way we combine elemental abundance 
limits to get constraints on cosmology and parti- 
cle physics. Partly for this reason, in the analysis 
to be described below the solar system D abun- 
dance of 2 X 10~^ is taken as a firm lower bound 
on D, and the previously quoted upper limit of 
10^^ is utilized as an upper limit on D+^He [43]. 
In the first place, it must be stressed that the 
Songaila et al result is still preliminary, and could 
easily be due to interloping hydrogen clouds along 
the line of sight of the system being observed. 
Also, other arguments based on galactic evolu- 
tion and the pre-solar D+^He abundance are in 
apparent contradiction with the result. In any 
case, if the Songaila et al result were to be con- 
firmed, it would require a discrete and dramatic 
shift in the entire BBN analysis, as we have de- 
scribed elsewhere. Thus, while one should keep 
in mind that such a parallel constraint space is 
a possibility, it is sufficiently "non-standard" so 
that it should be considered independent of the 
more generic systematic uncertainties we wish to 
concentrate on here. 

4.2. Results and Tables 

Tables 2-4 give our key results, adn the follow- 
ing description of them is taken from [3]. The 
data were obtained using 1000 Monte Carlo BBN 
runs at each value of rjiQ, with nuclear reaction 
rate input parameters chosen as Gaussian ran- 
dom variables with appropriate widths (see [1] 
for details) . In each case the number of runs 
which resulted in abundances which satisfied the 
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joint constraints obtained by using combinations 
of the upper limits on ^He, ^Li, and D+^He or the 
lower limit on D was determined. Limits on pa- 
rameters were determined by varying these until 
less than 50 runs out of 1000 (up to \//V statistical 
fluctuations) satisfied all of the constraints. 

Table 2 displays the upper limit on for var- 
ious values oiYp. As is shown, this was governed 
by the combination of "'He and D-|-'^He upper lim- 
its. Shown in the table are the number of accept- 
able runs out of 1000 when the two elemental 
bounds are considered separately and together, 
for an rjiQ range which was found to maximize 
the number of acceptable models. Throughout 
the "5 J,"'"'' region from .24 to .25, both the Yp and 
D-|-^He limits play a roughly equal role in deter- 
mining the maximum value of . We are able to 
find a remarkably good analytical fit for the max- 
imum value of Ni, as a function of Yp as follows: 

VV™"^ = 3.07 + 74:.Q7(Y^'"' - .240) (1) 

The linearity of this relation is striking over the 
whole region from .24 to .25 in spite of the in- 
terplay between the two different limits in deter- 
mining the constraint. Note also that this rela- 
tion differs from that quoted in [22] between Yp 
and in that the slope we find is about 13% 
less steep than that quoted there. The two for- 
mulae are not strictly equivalent in that the one 
presented in [22] presented the best fit value of 
Yp determined in terms of N^, while the present 
formula gives a relation between the maximum al- 
lowed values of these parameters, based on limits 
on the combination Yp and D+^He, and on the 
width of the predicted distribution. In this sense, 
eq. (1) is the appropriate relation to utilize when 
relating bounds on Yp to bounds on N^. 

Tables 3 and 4, which display the upper bounds 
on (]io- <ii'<^ perhaps even more enlightening. They 
demonstrate the sensitivity of the upper limit on 
rjiQ and hence flbaryon to the various other ele- 
mental upper limits as Yp is varied. Several fea- 
tures of the data are striking. First, note that 
^He completely dominates in the determination 
of the upper limit on (]io until Yp =.245, even for 
the most stringent chosen upper limit on ^Li. If 
this limit on ^Li is relaxed, then ''He dominates 
as long as the upper limit on Yp <.248! Also 



note that the "turn on" in significance of the Li 
contribution to the constrain is somewhat more 
gradual than the "turn off' of the ^He constraint. 
The former turns on over a range of rjiQ of about 
2, while the latter turns off over a range of about 
1-1.5. This gives one some idea of the size of the 
error introduced in determining upper bounds by 
using only either element alone, rather than the 
combination. Next, for a Yp upper limit which 
exceeds .248, the lower bound on D begins to be- 
come important. It quickly turns on in signifi- 
cance so that by the time the upper limit on Yp is 
increased to .25, ^He essentially no longer plays a 
role in bounding jjiq. Finally, note that both the 
relaxed bound on ^Li and the D bound converge 
in significance at about the same time, so that for 
rjiQ > 7.25, both constraints are significantly vio- 
lated. This implies a "safe" upper limit on rjiQ at 
this level, which corresponds to an upper bound 
^baryon < .163, assumiug a Hubble constant in 
excess of 40 km/sec/Mpc. We again stress that a 
value this large is only allowed if Yp exceeds .250. 
If, for example, Yp < .245, then the upper bound 
on ilhanjov is essentially completely determined 
by **He and is then at most 0.11. These limits 
may be compared to recent estimates of ^baryon 
based on X-ray determinations of the baryon frac- 
tion in clusters [49]. 

One final comment on the role of Yp in con- 
straining ?7io: It has been stressed that because 
of the logarithmic dependence of the former on 
the latter, that Yp cannot be effectively used to 
give a reliable upper bound on rjio. This is some- 
what deceptive, however. We can compare how 
much more sensitive the bound on jjiq is to Yp 
than the bound on N,j is by making a linear fit 
to the former relation and comparing it to (1). 
If we do this, we find first that the linear fit is 
quite good out to Yp as large as .245 (after which 
a quadratic fit remains good all the way out to 
.248, where the D and relaxed ^Li bounds begin 
to take over), and is given by 

jji"^"^ « 3.22 + 354(F™'"^ - .240) (2) 

Seen in these terms, the (]io upper limit is approx- 
imately 4.5 times more sensitive to the precise up- 
per limit chosen for Yp than is the upper limit. 
Thus, while there is no doubt that varying the up- 
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Table 2 

''Ho Abundanro Estimates <k N„ limits 



V 






# allowed 


models: 








{"He & [D+3He]}(-iHe:D+3He) 








mo=2.75 


2.80 


2.85 


2.90 


.240 


3.07 


40(603:148) 


52(429:254) 


46(268:376) 


38(170:534) 






??io=2-80 


2.85 


2.90 


2.95 


.241 


3.14 


38(532:171) 


46(354:309) 


39(219:470) 


35(131:625) 


.242 


3.21 


41(562:154) 


55(451:276) 


53(272:423) 


52(163:616) 


.243 


3.29 


17(588:110) 


32(410:220) 


46(266:378) 


36(184:513) 


.244 


3.36 


30(669:102) 


44(501:187) 


38(353:336) 


40(216:464) 






mo=2.85 


2.90 


2.95 


3.00 


.245 


3.43 


50(598:173) 


68(449:296) 


64(308:427) 


54(173:586) 


.247 


3.59 


27(635:84) 


30(480:184) 


47(338:306) 


39(185:488) 






??io=2-95 


3.00 


3.05 


3.10 


.250 


3.82 


45(491:207) 


47(364:374) 


50(225:495) 


32(131:587) 



Table 3 

''He and ' Li Abundance Estimates & jjiq limits 





(^Li_io < 2.3) 


# allowed models: 
{"He&^Li} (''He:'^Li) 


fLi.io < 5) 


# allowed models: 
{4He&^Li|(4He:^Li) 


.240 


3.26 


56 (60:998) 


3.26 


56 (60:1000) 


.241 


3.55 


45 (45:986) 


3.55 


45 (45:1000) 


.242 


3.89 


45 (47:905) 


3.89 


47 (47:1000) 


.243 


4.26 


50 (60:626) 


4.27 


46 (46:1000) 


.244 


4.64 


48 (92:296) 


4.71 


49 (49:1000) 


.245 


5.01 


45 (211:118) 


5.23 


62 (62:984) 


.246 


5.23 


51 (679:62) 


5.80 


46 (50:810) 


.247 


5.25 


52 (997:52) 


6.36 


48 (80:500) 



per limit on Yp has a more dramatic effect on the 
upper bound one might derive for rjiQ than it does 
for constraining Np,the quantitative nature of the 
relative sensitivities is perhaps displayed, for the 
relevant range of Yp, by comparing the linear ap- 
proximations presented here than by discussing 
logarithmic vs linear dependencies. More impor- 
tant, even recognizing the increased sensitivity of 
rjiQ on Yp, unless one is willing to accept the pos- 
sibility of a rigid upper bound on Yp greater than 
.247, it is overly conservative to ignore ''He when 
deriving BBN bounds on rjiQ. 

Finally, it is in this analysis that we rederived 
the minimum value of Yp such that BBN pre- 
dictions are consistent with observation. We ex- 
plored the range of 7710 allowed at the 95% con- 



fidence level (i.e. 50 out of 1000 models) as the 
value of Yp"''"' was reduced. For Yp < .239 no 
range of rjio was allowed when this constraint was 
combined with the D -|- ^He bound. Previously 
we derived a lower bound on 1',, of .238 if D+^He 
was used alone to first bound (]io. and then the 
rjiQ value was used to bound Yp (to compare to 
earlier such bounds (i.e. [11]). The new neu- 
tron half life would not change that bound. How- 
ever in any case the newly derived bound of .239 
obtained using the correlated constraints is more 
stringent, and more consistent. If the primordial 
helium abundance is determined empirically to 
be less than this value with great confidence, and 
the D + ^He upper limit remains stable, standard 
BBN would be inconsistent with observation. 
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Table 4 

^He, D and ^Li Estimates & rjio limits (^Li_io <5; -D-s > 2) 



y 

Prnax 




# allowed models: 






{"He & D & ^Li} ("He & D:"He & ^Li:D & 'Li) ("HeiDi^Li) 


.248 


6.94 


48 (136:53:156) (178:516:203) 


.249 


7.22 


52 (177:101:64) (654:217:136) 


.250 


7.24 


47 (191:113:47) (995:191:113) 



5. CONCLUSIONS AND CHALLENGES 

There can be little doubt that the present abil- 
ity of BBN to constrain cosmological parameters 
is almost completely governed by systematic un- 
certainties in our inferences of the actual light el- 
ement primordial abundances. Nevertheless, the 
fact that such systematic uncertainties need not 
be gaussian does not block our ability to utilize 
the statistically meaningful uncertainties in BBN 
predictions. As long as we are willing to cjuote 
conservative one-sided limits on the various abun- 
dances which incorporate reasonable estimates of 
the systematic uncertainties then the determina- 
tion of what confidence levels can be assigned to 
various theoretical predictions is straightforward. 
Moreover, as the observational limits on various 
elemental abundances is varied, the significance 
of the different elements for constraining cosmo- 
logical parameters varies. In addition, for a non- 
trivial range in jjiq, correlations exist between the 
various abundance predictions, and a self consis- 
tent use of all available constraints is important. 
Finally, Yp, in spite of its systematic uncertainty, 
plays a dominant role unless one is willing to ac- 
cept an upper limit of greater than .247. Beyond 
that, the convergence of D and ^Li limits suggest 
a safe upper bound of on the baryon density to- 
day of less than 16% of closure density. 

We thus find ourselves with on the verge of sev- 
eral possible interesting inconsistencies: 

1. If the primordial *Hc fraction can definitively 
be established to be less than 23.9% either D 
-|-^He estimates will have to be revised, or some 
more dramatic cosmological or particle physics- 
based alteration in BBN predictions will be re- 
quired 

2. Recent estimates of the baryon fraction in rich 
clusters (i.e. [49]) suggest that this fraction can 



be rather large. If the Universe is fiat, then the 
baryon fraction which one would derive based on 
the rich cluster estimates is at least a fraction of 
2 larger than that allowed by our BBN estimates, 
even with systematic uncertainties allowed for. 
Whether this is an indication that the universe 
is not fiat, or an indication that the cluster esti- 
mates themselves suffer from large possible sys- 
tematic uncertainties remains to be seen 

In any case, as time proceeds and more inde- 
pendent observations are made we will undoubt- 
edly get a better handle on the systematic un- 
certainties which presently limit the efficacy of 
BBN constraints. As I have shown here, several 
very interesting possible constraints on cosmol- 
ogy and particle physics are around the corner. 
Until then, BBN analyses still provide among 
the most useful tools to constrain fundamental 
physics. The updated tables and relations pre- 
sented here should provide a useful reference to 
allow researchers to translate their own limits 
on the light element abundances into meaningful 
bounds on and rjiQ. 

I would like to thank my present BBN collabora- 
tor Peter Kernan. for his significant contributions 
to the work discussed here 
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